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Abstract — In  the  study  of  mechanics  and  optimal  con¬ 
trol,  one  often  encounters  what  is  called  a  two-point 
boundary-value  problem  (TPBVP).  A  couple  of  meth¬ 
ods  exist  for  solving  these  problems,  such  as  the  Simple 
Shooting  Method  (SSM)  and  its  variation,  the  Multiple 
Shooting  Method  (MSM).  In  this  paper  a  new  method 
is  proposed  that  was  designed  from  the  favorable  as¬ 
pects  of  both  the  SSM  and  the  MSM.  The  Modified 
Simple  Shooting  Method  (MSSM)  sheds  undesirable 
aspects  of  both  previously  mentioned  methods  to  yield 
a  superior,  faster  method  for  solving  TPBVPs.  The 
convergence  of  the  MSSM  is  proven  under  mild  con¬ 
ditions  on  the  TPBVP.  A  comparison  of  the  MSM  and 
the  MSSM  is  made  for  a  problem  where  both  methods 
converge.  We  also  provide  a  second  example  where 
the  MSM  fails  to  converge  while  the  MSSM  converges 
rapidly. 


Although  not  the  first  to  investigate  the  solutions  of  TP¬ 
BVPs,  one  of  the  first  publications  to  approach  this  sub¬ 
ject  was  by  Keller  [2].  Those  initial  methods  were  and 
still  are  referred  to  as  shooting  methods. 

Keller  [3]  develops  the  SSM  and  the  MSM,  referring 
to  the  MSM  as  parallel  shooting,  and  also  proposes 
a  version  of  parallel  shooting  that  he  calls  "stabilized 
march.”  Several  years  later,  J.  Stoer  and  R.  Bulirsch  [5] 
explored  both  the  SSM  and  the  MSM  in  great  detail, 
while  providing  several  examples  and  hints  for  practi¬ 
cal  numerical  implementation. 

In  this  paper  a  new  method  is  proposed  for  the  solu¬ 
tion  of  two-point  boundary-value  problems  that  seems 
to  converge  faster  and  more  accurately  than  the  MSM. 
The  existence  and  uniqueness  of  solutions  to  the  TP¬ 
BVP  is  assumed. 
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1 .  Introduction 

A  general  TPBVP  can  be  written  in  the  following  form: 

y'(x)  =  f{x,y);  a<x<b  (1) 

r(y(a),y(b))  =  0,  (2) 

where  (2)  describes  the  boundary  conditions  satisfied 
by  the  system.  Examples  are  the  familiar  initial-value 
problem  (IVP)  and  first  order  necessary  conditions  ob¬ 
tained  by  an  application  of  the  Pontryagin  Maximum 
Principle  in  optimal  control  theory.  TPBVPs  from  op¬ 
timal  control  (unconstrained)  have  separated  boundary 
conditions  of  the  type  rx(y(a))  =  0  and  7*2  {y(b))  =  0. 
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Generally  speaking,  existence  and  uniqueness  theo¬ 
rems  for  two-point  boundary  value  problems  can  be 
quite  difficult;  however,  in  this  next  section  we  quote 
two  results  that  broach  this  topic.  The  first  is  from  Stoer 
and  Bulirsch  [5],  and  the  second  is  from  Keller  [2]. 

2.  Existence  and  Uniqueness 
Theorems 

An  existence  and  uniqueness  theorem  for  initial-value 
problems  can  be  found  in  Hale  [1],  which  is  but  one  of 
many  texts  that  provide  this  well  known  result.  On  the 
other  hand,  TPBVPs  may  have  multiple  or  no  solution 
at  all.  For  example,  consider  the  following  system: 


±l(t) 

1 

.  x2(t) 

g(x  1,2:2)  . 

where  g( ■)  is  a  continuous  function  of  its  arguments, 
and  xi(0)  ==  1,  xi(l)  =  —1.  One  can  easily  see  that 
the  first  equation  will  only  allow  values  of  x\  to  in¬ 
crease  as  time  increases.  Thus,  there  does  not  exist  a 
value  for  #2(0)  that  will  drive  the  value  of  x\  from  1 
at  t  —  0  to  —1  at  t  =  1.  Because  of  this  fact,  ex¬ 
istence  and  uniqueness  theory  for  TPBVPs  is  consid- 
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erably  less  developed  and  less  understood  than  that  of 
IVPs.  Despite  these  drawbacks,  below  are  two  exis¬ 
tence  and  uniqueness  theorems  that  are  applicable  to 
much  smaller  classes  of  functions  / ( x ,  y). 

General  Boundary  Conditions 

Theorem  2. 1:  For  the  two-point  boundary- value  prob¬ 
lem  (l)-{2),  let  the  following  assumptions  be  satisfied: 


2-  /fa*  y>  ?/)  satisfies  a  uniform  Lipschitz  condition 
on  R  in  y  and  yf. 


3.  /(e,|/,2/0  has  continuous  derivatives  on  D  which 
satisfy,  for  some  positive  constant  M, 


df 

(S)5T>0 


1.  /  and  Dyf  are  continuous  on  S  =  {(x,y)|a  <  x  < 
b,yefHn} 


(b) 


5/ 

dy' 


<  M 


2.  There  is  a  Jfc(  )  g  C[a,  b }  with  || Dvf(x,  y)||  <  k(x)  4‘  Tf,e  coefficients  in  (3)  satisfV  «o«i  >  0,  &0&1  > 

for  all  (x,  y)  g  5.  l^l  +  l^ol  7^  0 


3.  The  matrix 

P(u,v)  as  Dur(u,v) 

admits  for  all  u,  r;  €  3?n  a  representation  of  the  form 

P(ti,  v)  =  Po(7  +  Af  (u,  1/)) 

with  a  constant  nonsingular  matrix  P0  and  a  matrix 
M  =  M(u,  v),  and  there  are  constants  /i  and  m  with 

\\M(u,v)\\  <  p  <  1,  ||P0-1Dvr(u,r;)||  <  m 

for  all  €  9ftn. 

4.  There  is  a  number  A  >  0  with  A  +  /i  <  1  such  that 


Then  the  boundary  value  problem  (1)  has  exactly  one 
solution  y(x). 

Proof:  For  a  proof  of  this  theorem,  consult  Stoer 
and  Bulirsch  [5],  page  510.  H 

Separated  Boundary  Conditions 

For  a  theorem  that  will  apply  to  separated  boundary 
conditions,  we  consult  Keller  [2].  Consider  the  follow¬ 
ing  second-order  system: 

2/"  =  /(e,2/,J/');  a  <  x  <  b 

aoy{a)  -  aiy'(a)  =  a,  \a0\  +  \ax\  £  0; 

boy(b)  +  &ij/'(6)  =  /?,  |60|  +  |&i|  0.  (3) 


Then  the  boundary-value  problem  in  (3)  has  a  unique 
solution. 

Proof:  For  proof  of  this  theorem,  consult  Keller 
[2],  page  9.  ■ 

Remark.  Assumptions  1-4  of  Theorem  2.1  are  very 
restrictive  sufficient  conditions.  Even  simple  boundary 
conditions  exist  that  do  not  satisfy  assumption  3;  such 
is  the  case  with  separated  boundary  conditions. 

Optimal  Control 

Now  consider  the  optimal  control  problem  of  finding  a 
ti(-)  for  the  following  system: 

x  =  /(i,u),  x(0)  =  I0)  a:(l)  =X1,  (4) 

such  that 

J{u)~  f  L(x,u)cff 
Jo 

is  minimized.  The  Pontryagin  Maximum  Principle 
yields  the  existence  of  functions 

P(t)  =  \Pl(t)p2(t)  p„(*)]T 

witht  €  [0,1];  H(x,u,p)  =  L(x,u)  +pTf(x,  it),  and 
u*  =  arg  minu  H(x,  u,  p)  such  that 

[?](*)-[-/ 

satisfies  x(0)  =  x0  and  x(l)  =  xx .  If 

VH(x)o,rgm\nuH(x)u1p)ip)  is  Lipschitz  continu¬ 
ous  in  the  x  and  p  variables  then  we  have  uniqueness. 
A  sufficient  condition  is  the  twice  differentiability  of 
/f(x,u,p). 


Theorem  2.2:  Let  the  function  f(x ,  y,  y')  in  (3)  satisfy 
the  all  of  the  following: 

1 .  /(x,  y,  y')  is  continuous  on  D  ~  {(x,  y,  y')  | 
a  <  x  <  6,  y2  +  (y')2  <  oo} 


Now  that  we  have  proof  of  existence  and  uniqueness 
for  small  classes  of  TPBVPs,  let’s  explore  the  current 
methods  commonly  used  to  numerically  solve  such  a 
problem  and  take  a  look  at  the  new  method  that  we 
propose. 
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3.  Current  Methods 

Although  Theorem  2.1  does  not  apply  to  the  case  of 
separated  boundary  conditions  and  Theorem  2.2  itself 
may  be  somewhat  restrictive,  separated  boundary  con¬ 
ditions  are  the  most  commonly  encountered  in  optimal 
control.  Because  of  this,  separated  boundary  condi¬ 
tions  will  be  used  for  explanation  purposes.  The  system 
now  becomes 

y'(x)  =  f(x,y );  a  <  x  <  b 
Ay(a)  =  a ,  By(b)  =  0,  (5) 

where  A  and  B  sltg  m  x  n  matrices  with  rank(A)  + 
rank(B)  =  n. 

Simple  Shooting 

The  Simple  Shooting  Method,  as  the  name  implies, 
is  the  simplest  method  of  finding  a  solution  to  such  a 
problem.  The  idea  is  to  convert  (5)  into  an  initial-value 
problem  (IVP): 

y'(x)  =  f(x,y)]  a  <  x  <b 

y(a)  =  Va,  (6) 

where  ya  is  composed  of  known  states  from  Ay  (a)  = 
a  and  guesses  for  the  unknown  states  s0.  Now,  y(x)  6 
ol  6  and/?  €  $ftp.  A  necessary  condition 
to  keep  the  problem  from  being  inconsistent  is  that 
m  +  p  =  n.  To  form  an  IVP  out  of  (5),  one  needs 
to  guess  initial  conditions  for  the  (n  —  m)  components 
of  y(x)  that  do  not  already  have  initial  conditions  at 
x  —  a.  Let  so  G  5ftn“m  be  the  guess  for  the  unknown 
initial  conditions  and  s*;  k  >  1  subsequent  corrections 
of  the  vector  sq.  With  sq,  one  can  now  integrate  (6)  for¬ 
ward  in  the  time  variable  x.  Integration  is  performed 
from  x  —  a  all  the  way  to  x  =  b.  Then  compute  the 
error  e  =  ||i?2/(b)  -  /?||2.  With  this  information,  a  cor¬ 
rection  is  made  to  the  initial  guess  s0  to  yield  su  and 
the  integration  is  performed  again.  This  process  is  re¬ 
peated  over  and  over  until  e  <e,  where  e  >  0  is  small. 
How  the  correction  is  made  will  be  addressed  shortly. 

For  illustration  purposes  only,  consider  Figure  1  which 
represents  an  example  of  Simple  Shooting  for  (6).  We 
assume  that  there  exists  a  unique  solution  to  the  prob¬ 
lem.  Every  point  on  the  plot  represents  a  vector  in 

There  can  be  serious  problems  with  the  accuracy  of  the 
SSM.  The  problems  occur  when  making  the  correction 
to  the  sjk  vector.  This  vector  is  usually  corrected  us¬ 
ing  a  modified  Newton’s  Method,  and  in  practice  the 
system  must  be  linearized  to  use  this  method.  If  e  is 
large,  then  convergence  can  be  quite  slow  (please  re- 


Figure  1.  Illustration  of  the  Simple  Shooting  Method 


fer  to  page  51 1  of  Stoer  and  Bulirsch  [5]).  This  draw¬ 
back  of  the  SSM  can  be  fixed  by  implementing  what  is 
known  as  the  Multiple  Shooting  Method. 

Multiple  Shooting 

The  Multiple  Shooting  Method  begins  with  the  choice 
of  a  Lipschitz  continuous  function  <p(x)  that  satisfies 
A<p(x )  =  a  and  B<p(x)  —  f3.  An  initial  guess  of  un¬ 
knowns,  so,  must  be  made.  Then,  (6)  is  integrated  until 
II y(x>  so)  -  y>(z)||  >  £  for  some  e  >  0.  We  designate 
the  time  variable  at  this  point  as  x\ .  Now  the  integra¬ 
tion  of  the  system  continues  from  x±  using  p(x i)  as 
the  initial  ’guess’  for  the  solution.  This  process  contin¬ 
ues  until  the  integration  reaches  x  =  b.  Now  the  error 
function  e(s)  =  §y(xi)  -  2/i||2  is  formulated  where 
s  =  [so  2/i  •  •  •  yk-i]T  and  y{  is  the  initial  state  for 
the  trajectory  in  the  interval  [a^Xi+i].  After  this  is 
accomplished,  a  correction  is  made  to  s  using  a  modi¬ 
fied  Newton’s  method,  and  the  process  is  repeated.  The 
starting  trajectory  is  not  used  after  the  first  iteration. 

Figure  2  illustrates  the  MSM.  Once  again,  three  iter¬ 
ations  are  displayed  for  this  method.  The  correction 
process  stops  when  e(s)  <  ei  <  e  for  some  ex  >  0. 

A  problem  with  the  MSM  is  the  discontinuity  of  the 
trajectory  found  by  the  MSM  at  the  points  Xjj  i  = 
lj  •  •  •  >  &  “  1.  The  integration  and  corrections  of  s 
will  continue  until  a  desired  level  of  closeness  is  de¬ 
termined,  but  this  final  value  of  the  vector  s  can  still 
be  far  from  an  optimal  solution  due  to  the  unstable  na¬ 
ture  of  many  systems  in  the  forward  direction.  If  (6) 
is  re-integrated  to  result  in  one  continuous  trajectory 
for  the  system,  the  end  values  By(b)  need  not  be  any¬ 
where  close  to  p  and  almost  certainly  will  not  be.  An¬ 
other  problem  is  computation.  During  the  process,  one 
must  invert  many  matrices  of  the  size  [n(k  -  1)  +  p], 


3 


where  k  can  be  quite  large  depending  on  the  guesses 
[so  <p{x  1)  •  •  •  ip{xk-i)]T  for  the  initial  trajectory. 
Note  that  k  cannot  be  reduced  even  as  the  guesses  im¬ 
prove. 

Example  3.1 :  Consider  the  following  system: 


where  0  <  x  <  1.  This  system  was  solved  with  the 
‘bad*  initial  guess  s0  =  [-100  2]T  with  the  parame¬ 
ters  of  the  code  as  follows: 

•  The  time  step  h  —  0.01. 

•  e  was  set  to  1. 

•  Newton’s  Method  is  stopped  when  the  solution  is 
found  to  be  in  an  ex  ball  of  size  ex  —  10"3 

For  this  example  (p(x )  =  x[l  1]T  +  [1  1]T.  The 
following  results  were  obtained.  The  MSM  corrects 
s0  to  [0.70  0.70]T  in  8.543  seconds.  Figure  3  shows 
the  discontinuous  segments  obtained  after  the  conver¬ 
gence  of  the  iterations.  One  can  see  the  discontinuous 
segments  connecting  the  given  initial  and  desired  final 
states.  However,  when  this  very  system  is  re-integrated 
using  the  $o  vector  determined  to  be  the  correct  one  by 
the  MSM,  the  first  two  components  of  the  solution  do 
not  end  up  within  10“ 3  of  [2  2]T  (please  see  Figure  4. 
Instead,  the  solution  ends  up  being  close  to  [2.5  2.5]T. 
Obviously,  these  are  not  desirable  results. 

Note  that  while  using  the  MSM  one  does  not  have  con¬ 


trol  over  the  error  in  the  states  at  the  final  time.  It  de¬ 
pends  on  the  particular  system  being  considered.  For 
the  sample  problem,  in  order  for  the  error  between  the 
actual  and  desired  final  states  to  be  small,  we  reduced 
the  time  step  to  h  —  0.001.  This  time  MSM  corrects 
the  unknown  states  so  to  [0.42  0.42]T  in  16.634  sec¬ 
onds.  Figure  5  shows  both  the  result  of  MSM  and  the 
final  trajectory  after  re-integration. 


Figure  3.  Discontinuous  segments  connecting  the  ini¬ 
tial  and  desired  final  states  while  using  the  MSM. 


Figure  4.  Plot  of  the  states  while  using  the  result  of 
the  MSM  ( h  =  0.01  seconds). 


4.  Modified  Simple  Shooting 

Description  of  Algorithm 

The  Modified  Simple  Shooting  Method  begins  with  the 
selection  of  a  Lipschitz  continuous  starting  path  of  in¬ 
tegration,  <p(x) ,  such  that  Aip(a)  —  a  and  B<p(b)  —  (3. 
Again,  an  initial  guess  of  unknowns,  so  must  be  made. 
Then,  (6)  is  integrated  until  ||y(x,  s0)  -  y?(:r)||  >  £  f°r 
some  e  >  0.  Then  we  designate  sXo  =  sq.  A  modified 
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Figure  5.  Plot  of  the  states  while  using  the  result  of 
the  MSM  ( h  =  0.001  seconds). 


rigure  j  musiraies  me  Moamea  simple  snooting 
Method.  In  this  case,  it  took  three  overall  ’shots’  to 
integrate  from  x  =  a  to  x  ~  b. 


Figure  6.  Illustration  of  the  Modified  Simple  Shooting 
Method 


Newton’s  Method  is  then  used  to  correct  the  ’guess’ 
slo.  The  iteration  stops  when  || y(x>  sik)  —  <p(x) ||  <  £\ 
where  ex  is  chosen  such  that  £\  <  e.  We  then  let 
s\  =  sik ,  and  proceed  with  the  integration  of  the 
system,(6)  where  y{a)  is  found  using  Ay(a)  —  a  and 
«i. 

The  modified  Newton’s  Method  mentioned  above  is 
found  in  Stoer  and  Bulirsch  [5].  The  following  is  an 
outline  of  that  method  for  the  first  iteration: 

1.  Choose  a  starting  point  s0  € 

2.  For  each  i  —  0, 1, . . .  define  $ii+1  from  Si.  as  fol¬ 
lows: 

(a)  Set 

di  =  DF(sli)-1F(sli), 

1 

7i  “  cond(I>F(si4))  ’ 

and  let  hi(r)  =  h(sxi  —  rd*),  where  h(s)  = 
F(s)tF(s).  Determine  the  smallest  integer  j  >  0  sat¬ 
isfying 

/li(2-J)</li(0)-2^|||di||  \\Dh(Sli)\\. 

(b)  Determine  so  that  h(sxi+l )  =  mino<K<j  h{( 2-K), 
and  let 

Sli+-1  =  5ii  +  \{di. 

The  MSSM  continues  until  ||2/(x,sg)  -  y?(x)||  <  e  at 
each  x  G  [a,  b].  In  this  last  step,  we  are  performing  ex¬ 
actly  the  SSM  for  the  original  system,  but  with  a  start¬ 
ing  initial  guess  sq  that  keeps  By  (b)  close  to  ft.  This 
will  prevent  any  numerical  divergence. 


Convergence  of  the  MSSM 

Theorem  2.1  provided  us  with  a  somewhat  limited  exis¬ 
tence  and  uniqueness  theorem.  Theorem  2.2  was  more 
useful  in  the  sense  that  it  applied  to  the  case  of  sep¬ 
arated  boundary  conditions.  For  the  purpose  of  this 
section,  we  shall  assume  the  existence  of  a  solution  to 
(5)  for  an  entire  family  of  boundary  conditions.  More 
specifically,  suppose  that  the  BVPs 

y'{*)  = 

Ay{a)  —  a,  By(x)  =  Bip(x) 

have  unique  solutions  where  </?(*)  was  the  function  cho¬ 
sen  initially.  This  is  necessary  to  assume,  as  a  solution 
to  the  overall  problem  may  not  directly  imply  the  exis¬ 
tence  of  a  solution  to  one  of  the  intermediate  reduced 
problems.  Now  we  can  consider  the  main  theorem  of 
this  report  on  the  convergence  of  the  Modified  Simple 
Shooting  Method  with  certain  assumptions. 

Theorem  4.1:  Consider  the  Two-Point  Boundary- Value 
Problem  as  described  in  (5).  Let  y(x)  denote  the  solu¬ 
tion  to  this  problem.  The  Modified  Simple  Shooting 
Method,  as  described  earlier,  converges  to  y(x)  when 
applied  to  (5). 

Proof:  In  order  for  the  MSSM  to  converge  to 
y(x),  it  must  first  converge  to  yi(x)  at  each  interme¬ 
diate  point  X{ ,  where  i  =  1, 2, . . . ,  k  —  1,  and  yfx) 
is  the  solution  of  the  reduced  problem  on  [a,  X{]  for 
i  =  1,2,  ...,&-l.  For  example,  yi(x)  is  the  solution 
to  the  problem 

y'lix)  =  f{x,yi),  a<x<  Xi 
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Ay(a)  =  q,  By(xi)  =  <p(x i) 


where  p(x)  is  the  reference  path  mentioned  in  the  de¬ 
scription  of  the  algorithm.  As  such,  it  is  only  necessary 
to  show  two  things  to  complete  this  proof. 

1.  The  sequence  of  points  {xn}™=l  e  9?  converges  to 
the  right  endpoint  b. 

2.  The  SSM  converges  when  existence  is  known. 

The  latter  is  a  result  of  the  modified  Newton’s  Method, 
which  is  guaranteed  to  find  a  solution  for  large  classes 
of  functions,  if  it  exists.  As  existence  of  a  solution  y  (x) 
is  being  assumed,  the  modified  Newton’s  Method  guar¬ 
antees  that  the  SSM  converges  to  y(x). 

Now  assume  that  after  the  modified  Newton’s  Method 
is  performed  at  x  €  [a,  6],  the  integration  proceeds 
to  x*  >  x.  It  is  necessary  to  show  that  (x*  -  x) 
is  bounded  below  by  a  positive  number.  y(x,sm)  is 
a  continuous  function  of  x  and  p(x)  is  a  Lipschitz 
continuous  function  of  x.  Thus,  d((p;y,sm)(x)  = 
l|y(x>Sm)  —  ¥>(x)||  is  a  Lipschitz  continuous  func¬ 
tion  of  x.  By  the  compactness  of  [a,  b],  there  is  a 
uniform  Lipschitz  constant,  k  for  d(<p;  y,  •).  Thus, 
\d(<p;y1sm)(x)~d((p]y,sm)(x,*l)\  <  fc|x-x*|.  When 
the  integration  of  the  system  stops,  sm)(x)  - 

d(<p\y,sm)(x*)\  =  e— ei.Thene— ei  <  fc(x*-x), or 
(x*  —  x)  >  >  0.  This  means  that  the  integration 

continues  past  x,  and  by  the  compactness  of  [a,  6],  the 
process  must  extend  all  the  way  to  x  =  b.  ■ 

Examples 

Remark.  All  computations  in  the  following  examples 
were  performed  in  the  MATLAB  environment.  Ver¬ 
sion  6.1.0.450  Release  12.1,  running  on  a  Microsoft 
Windows  2000  Professional  operating  system  with  an 
AMD  Athlon  processor  running  at  1.2  GHz. 

A  Linear  Example- 

Example  4.1:  Consider  the  simple  system  of  Example 
3.1.  The  MSSM  was  applied  to  the  same  system  with 
the  same  parameters  in  the  code. 

The  MSSM  finds  a  solution  in  1.001  seconds  with  s 
corrected  to  [0.39  0.39]r,  which  is  about  8  times  faster 
than  the  MSM.  The  corrections  to  s  by  both  methods 
do  not  give  the  same  result.  The  discontinuity  problem 
does  not  arise  with  the  MSSM,  since  the  final  trajectory 
is  a  continuous  one. 

An  Application  in  50(3)— Now  we  will  look  at  an 
example  in  the  three  dimensional  special  orthogonal 


Figure  7.  Plot  of  the  states  for  the  trajectory  planning 
problem  using  MSSM. 


group  or  50(3).  50 (n)  is  defined  as  follows: 


SO(n)  =  {A  €  XnXn\detA  =  1,  AT A  =  /nxn} . 

This  example  is  actually  an  optimization  problem  on 
the  set  of  orientation  matrices  in  three  dimensional 
space.  Consider  the  following  system. 

Q  =  QCl 


'  n  ' 

Pi 

A 

= 

-5P1  xfi  +  fix  i(P2  X  ft) 

.A . 

-5P2  X  ft-Pj 

Q  j  Piy  and  P2  are  all  functions  of  time,  but  the  de¬ 
pendence  on  time,  £,  has  been  suppressed  for  ease  of 
writing.  ft,  Pi,  and  P2  are  vectors  in  9E3,  whereas  Q  is 
a  matrix  in  9ft3*3.  Cl  is  a  little  more  complicated;  it  is 
a  skew-symmetric  matrix  formed  from  the  vector  ft  as 
such. 

If 


ft  = 


ftj 

^2 

ft3 


then 


ft  = 


0  ”“ft3  ft2 

ft3  0  — -fti 

— ft2  ftl  0 


Qinitiah  ^initial*  Q desired  and  ft^eatred  are  known. 
The  solution  is  sought  for  0  <  t  <  1. 


For  this  particular  example,  further  obstacles  were  to 
be  overcome.  It  was  required  that  Q(t)  e  50(3)  at  ev¬ 
ery  instant  of  time  t .  Furthermore,  there  was  the  obsta¬ 
cle  that  we  must  integrate  forward  in  time  four  equa¬ 
tions,  three  of  which  are  vector  valued  equations  and 
one  matrix  valued  equation.  But,  the  matrix  equation 
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depended  only  upon  the  value  of  ft  and  hence  Q  at  each 
instant  of  time.  Because  of  this  and  the  fact  that  none 
of  the  other  equations  depend  on  <2,  it  was  possible  to 
integrate  the  matrix  equation  separately,  but  still  for¬ 
ward  in  time  at  the  same  time  as  the  vector  equations. 
To  keep  Q(t )  €  50(3),  Rodriques’  formula  was  used, 
which  can  be  found  in  Murray,  Li,  and  Sastry  [4].  The 
formula  is^ 

ene  —  /  +  Q  sin  e  +  ft2(i  —  cosO). 

The  initial  and  final  values  of  Q  and  ft  were  generated 
randomly  by  MATLAB  and  were  exactly  the  same  for 
both  the  MSSM  and  the  MSM.  Those  values  were 

0.57  0.82  0.08 

Qimtiai=  -0.78  0.57  -0.25  , 

-0.25  0.08  0.97 

'  0.95  ‘ 

ft initial  —  4.34  , 

7.09  _ 

’  -0.27  0.96  -0.01  ' 

Q desired  ~  “0.84  -0.23  0.50  , 

0.48  0.14  0.87 

'  1.90  " 

ft desired  ~  8.67  . 

_  4.18  _ 

The  results  for  this  example  again  heavily  favored  the 
Modified  Simple  Shooting  Method.  After  47.12  sec¬ 
onds,  the  MSSM  obtained  the  following  values  for  fi 
and  Q  at  time  t~  1, 

|  -0.27  0.96  -0.01  ' 

Q(l)  =  -0.84  -0.23  0.50  , 

0.48  0.14  0.87 

’  1.90  ’ 

ft(l)  =  8.67  . 

4  19  . 

The  Multiple  Shooting  Method  proved  inadequate  for 
this  problem.  After  606.89  seconds,  the  MSM  obtained 
these  results, 

'  0.09  -0.01  1.00 

Q(l)  =  0.99  -0.06  -0.09  , 

0.06  1.00  0.00 

2.10 

0(1)  =  10.81  . 

5.89 

To  measure  closeness  the  normal  Euclidean  norm  can 
be  used  to  compare  ^desired  and  ft(l);  however,  this 


is  not  the  case  with  Qdesired  and  Q(l).  One  must  be 
more  careful.  It  is  desired  to  have  a  measure  of  close¬ 
ness  within  the  group  50(3),  not  the  space  of  all  3  x  3 
matrices.  To  do  so,  we  take  the  matrix  logarithm  of 
the  quantity  Q(l)T Qdesired,  which  will  yield  a  skew- 
symmetric  matrix.  We  then  take  the  norm  of  this  ma¬ 
trix.  The  Multiple  Shooting  Method  yields 

|| ^desired  ~  ft(l)||  =  2.74, 

II  log(Q(l)TQdestred)||  —  2.62. 

The  Modified  Simple  Shooting  Method  yields 

IliWei  -  n(i)||  =  0.0013, 

II  log(Q(l)TQdesired)||  =  0.0003. 

The  MSM  took  more  than  10  times  as  long  and  the 
results  are  not  close  at  all  to  the  desired  values;  these 
results  speak  for  themselves. 
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Figure  8.  Plot  of  ZYX  Euler  angles. 
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Figure  9.  Q  matrix  components. 
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Figure  10.  Angular  velocity. 


5.  Conclusions 

In  this  report,  a  new  method  for  solving  two-point 
boundary-value  problems  was  described.  Although 
convergence  is  not  so  laborious  to  investigate,  it  was 
shown  by  two  examples  of  theorems  how  difficult  it 
is  to  prove  existence  and  uniqueness  for  two-point 
boundary- value  problems.  Nonetheless,  three  exam¬ 
ples  were  given,  among  many  performed,  that  clearly 
show  that  the  Modified  Simple  Shooting  Method  per¬ 
forms  better  and  faster  than  the  Multiple  Shooting 
Method. 

First,  it  requires  the  inversion  of  much  smaller  matri¬ 
ces  than  those  required  to  be  inverted  in  the  MSM.  The 
MSSM  requires  the  inversion  of  matrices  that  are  n  x  n; 
whereas,  the  MSM  requires  the  inversion  of  matrices 
that  are  [n(fc  - 1)  -f  p]  x  [n(fc  - 1)  +p].  This  fact  alone 
could  account  for  many  seconds  of  computation  time 
saved  as  systems  become  larger  and  larger. 

Another  fact  that  makes  the  MSSM  more  appealing 
is  continuity  of  integration  trajectory.  This  property 
is  very  important  in  optimal  control  problems  where 
the  systems  are  unstable  in  forward  tipie.  The  Modi¬ 
fied  Simple  Shooting  Method  integrates  the  system  in 
one  continuous  path  every  time  it  shoots  the  system  for 
an  updated  s  vector.  The  Multiple  Shooting  Method 
does  not  have  this  characteristic.  In  fact,  for  a  partic¬ 
ular  example  of  the  MSM,  if  k  intermediate  shots  are 
taken  then  every  overall  shot  of  the  system  from  x  =  a 
to  x  —  b  will  consist  of  k  -  1  discontinuities.  Each 
of  these  discontinuities  are  impossible  to  correct.  The 
best  the  method  has  to  offer  is  to  reduce  the  magni¬ 
tude  of  the  discontinuities.  Due  to  the  instability  of 
many  systems  in  the  forward  direction,  an  erroneous 
solution  may  be  obtained  if  the  system  is  re-integrated 


with  what  is  supposed  to  be  an  accurate  approximation 
to  the  actual  unknown  initial  conditions. 
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